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ABSTRACT 

We have used the integral field spectrograph OSIRIS and laser guide star adaptive optics at Keck 
Observatory to obtain high angular resolution (0'.'06), moderate spectral resolution (R ~ 3800) images 
of the bipolar jet from the Herbig Ae star LkHa 233, seen in near-IR [Fe II] emission at 1.600 & 1.644 
/im. This jet is narrow and tightly collimated, with an opening angle of only 9 degrees, and has an 
average radial velocity of ~ 100 km s . The jet and counterjet are asymmetric, with the red-shifted 
jet much clumpier than its counterpart at the angular resolution of our observations. The observed 
properties are in general similar to jets seen around T Tauri stars, though it has a relatively large 
mass flux of 1.2 ± 0.3 x 1O _7 M0 year -1 , near the high end of the observed mass flux range around T 
Tauri stars. We also spatially resolve an inclined circumstellar disk around LkHa 233, which obscures 
the star from direct view. By comparison with numerical radiative transfer disk models, we estimate 
the disk midplane to be inclined i = 65 ± 5° relative to the plane of the sky. Since the star is seen only 
in scattered light at near- infrared wavelengths, we detect only a small fraction of its intrinsic flux. 
Because previous estimates of its stellar properties did not account for this, either LkHa 233 must 
be located closer than the previously believed, or its true luminosity must be greater than previously 
supposed, consistent with its being a ~ 4M Q star near the stellar birthline. 

Subject headings: ISM: jets and outflows — ISM: Herbig-Haro objects — stars: individual (LkHa 
233) — stars:pre-main-sequence 
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1. INTRODUCTION 

Bipolar outflows from young stars feature prominently 
in many of the most spect acular images of our universe 
(see, e.g. iBallv et al]|2007|) . But far from being merely 
aesthetically pleasing, these jets play crucial roles in star 
formation. Molecular cloud material must shed most of 
its angular momentum before it can accrete onto a new- 
born star, and jets have been identified as a key mech- 
anism for this, removing angular momentum from disks 
and allowing accretion to continue (|Rav et al.ll2007l and 
references therein) . Recent observations have for the first 
time p rovided indications of jet rotation ( Bacciott i et al 



time provided indications or ict rotation (Bacciotti ct al 
2002; Cof fev et alJl2004t iWoitas et al.ll2005H Coffcv et a' 



20071) . suggesting that outflows may indeed carry angu- 
lar momentum away from their origin. Jets have also 
been implicated in regulating the overall efficiency of star 
formation by injecting tu rbulence into molecular clouds 
dMatzner fe McKedl2000h . 

The physical processes which drive these outflows re- 
main poorly understood. Several competing theories 
have been proposed to explain how gas is accelerated 
and collimated. The leading contenders all invoke mag- 
netic fo rces to shape the outflows, with the "X-wind" 
model (jShang et al.l l2007t ) positing acceleration occurs 
near the star at the radius of mag netospherical trunca- 
tion, while the "disk wind" model (jPudritz et al.ll2007h 

Electronic address: mperrin@astro.berkeley.edu 

1 Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 



posits acceleration from a wide range of radii across the 
circumstellar disk. Alternatively, outflows could origi- 
nate in spherical coronae, whic h are then collim ated by 
surrounding magnetic fields (jSautv et al.l Il999t ) . Out- 
flows may even evolve th rough each of these processes in 
turn (|Sautv et al.| [2003). Distinguishing between these 
theories will require detailed studies of outflows as close 
to their origins as possible. At present, high angular res- 
olution observations of outflows have mostly focused on 
older T Tauri s tars approaching the main sequence (see 
iRav et al"1l2007L and references therein). Observations of 
more sources at younger ages are needed, especially con- 
sidering that outflows are intimately rel ated to accretion, 
which is greatest at the youngest ages (jHartmann et al.l 
fl99l . 

In addition to age, we also need to understand how 
outflows vary with stellar mass. Because jets are thought 
to be driven by magnetic fields, we might suppose that 
their properties should differ between fully convective T 
Tauri stars and more massive stars which lack a surface 
convection zone. Indeed, outflows from high mass YSOs 
are much less collimated than jets f rom T Tauri stars , 
with opening angles of 30 — 60° (e.g. iHunter et al.lll997l : 
IShepherd etafl 1 19981 I2Q01). The Herb ig Ae/Be stars 
Verbid ll960TlWaters fc Waelkensl fl998h are intermedi- 
ate between these two regimes, so studies of their out- 
flows can potentially clarify this transition in jet prop- 
erties, which in turn may shed light on the underlying 
physics through which jets are produced. Yet relatively 
few outflows from Herbig Ae/Be stars ha ve been studied 
in de tail, particularly at high resolution (jMundt fe Ravi 
fl99l . In a very few cases, studies of parsec-scale 
outflows from Herbig Ae stars have suggested that on 




Perrin and Graham 




o 

Arcsec 



Fig. 1. — Left: This J,H,K S composite from the OSIRIS imager shows the X-shaped envelope surrounding LkHc* 233. The dashed 
rectangle shows the field zoomed in to on the right, while the solid line indicates the outflow axis. Right: A portion of the OSIRIS IFU 
data cube, centered in wavelength on the [Fe II] 1.644 fim line (in green) and rotated to show the jet axis vertical. This image reveals a 
tightly collimated bipolar outflow from LkHa 233 perpendicular to an edge-on circumstellar disk. The jet is undetectable in the broadband 
image; only with the spectral resolution of OSIRIS can it be detected. In the right hand panel, the blue-shifted jet points down and the 
red-shifted jet is upward. In these data we also see scattered light from a roughly edge-on circumstellar disk and a larger bipolar nebula 
extending perpendicular to the disk pla ne. The jet has excavated a cavity i n this envelope, the limb-brightened edges of which form the 
well-known "X" centered on LkHa 233 (Aspin et al. 1985; Perrin ct al. 2004). 



large spatial scales their velocities and collimation are 
largely consistent with those observed in T Tauri jets 
(jMcGroartv et al.ll2004h . 

But a lack of high angular resolution studies of Her- 
big Ae star outflows has thus far prevented detailed 
comparison with outflows from T Tauri stars. For T 
Tauri stars, we now have a sophisticated understand- 
ing of outflow properties on scales of a few tens of AU, 
derived from high angular resolution observations. Be- 
cause stellar outflows are detectable primarily in emis- 
sion lines from shock-cxcitcd ions, line ratios can provide 
measurements of many physical parameters within the 
outflow. When both spatially and spectrally resolved 



data are available, ven 
have been obtained (e.g. 


i detailed portraits of outflows 


Bacciotti et al.lll999tlBeck et all 


12004 lHartigan & Morse 


20071). Since outflows have char- 



acteristic spatial scales of a few AU (corresponding to 
O'.'Ol — O'.'l at typical distances) and velocities of tens to 
hundreds of km s , the necessary observations combine 
high angular resolution (< O'.'l) with moderate spectral 
resolution (R ~ 2000 - 4000). Having both spatially 
and spectrally resolved data is key to disentangling the 
outflow structures as close to the disk as possible. In 
optical wavelengths, such as for the [S II] 6717, 6731 
A and [O I] 6300, 6364 A emission lines, the neces- 
sary angular re solution is available o nly frorn space (e.g.. 



Hartigan et al.lll999t iBacciotti et aTll2000t iWoitas et al l 
2002t ICoffev et al.ll2004 lHartigan fe Morsel I2007D How- 



ever, the arrival of integral field spectroscopy with adap- 
tive optics on 8-10 m telescopes has now opened the 
door to high angular resolution studies of jets from the 
ground. At the near-infrared wavelengths suitable for 
adaptive optics, several emission lines of [Fe II] be- 
tween 1.2 — 1.6^m ar e ideal tracers of shocked outflows 
(jGraham et alj |1987[ ). Furthermore, these lines pro- 



vide diagnostic potential enabling t he determination of 
physical conditions w ithin the jet (|Pesenti et al.l 120031 : 
lHartigan et al.ll2004h . and ultimately the mass flux, a 
key parameter in any model of outflow physics. 

In this paper, we seek to investigate whether the sim- 
ilarity on large spatial scales between outflows from T 
Tauri and Herbig Ae stars still holds true on finer spa- 
tial scales. Toward that end, we present the first high 
angular resolution integral field spectroscopy of an out- 
flow from a Herbig Ae star, enabling us to investigate its 
collimation, kinematics, and certain physical parameters 
on angular scales of ~ 60 milliarcseconds. By compari- 
son with the properties of jets seen from T Tauri stars, 
we hope to clarify the influence of stellar mass upon out- 
flow physics. Also, from a technical perspective, we wish 
to demonstrate the utilty of laser guide star adaptive 
optics (LGS AO) for obtaining high angular resolution 
spectroscopy of an outflow from a fairly faint and ob- 
scured source. Because many young stars which drive 
jets are still deeply embedded in their dusty birthplaces, 
only a limited number are accessible to traditional natu- 
ral guide star AO. The growing availability of LGS AO 
on large telescopes promises to enable high angular res- 
olution studies of a greatly increased target sample, par- 
ticularly at the youngest ages. 

We present here observations made with Keek's new in- 
tegral field spectrograph OSIRIS, using laser guide star 
adaptive optics, which show in exquisite detail a tightly 
collimated jet from the Herbig Ae star LkHa 233 (also 
V375 Lac, HBC 313). T his star, part of the original Ae 
star sample identified by (|HerbigHl960T ). is one of the rel- 
atively few Herbig Ae/Be stars known to possess a jet. 
It is an A4 star, with logT e // = 3.93 (| Hernandez et all 
l2004h . Its distance has genera lly been considered to be 
880 pc (jCalvet fe Cohenlll978D . We adopt this distance 
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Fig. 2. — One dimensional spectra extracted from the OSIRIS data cube. Left: Square apertures indicate the regions extracted to 
produce the spectra at right, shown in corresponding colors. The red aperture is centered on knot B. Center and Right: One dimensional 
spectra extracted from the data cube for those apertures, in OSIRIS' Hn3 and Kn3 filters; the light grey shading shows the statistical 
uncertainty. For clarity, the emission line spectra have been scaled up xlO and x20 relative to the continuum. The outflow is clearly 
detected in both the 1.600 and 1.644 (im lines of [Fe II] , and excess emission in Br7 is seen at the location of the star. There is a hint 
of 2.122 (im H2 emission, but it unfortunately falls at the very edge of our bandpass so its reality is uncertain. The four features marked 
"X" are spectral ghosts caused by a grating misalignment during OSIRIS commissioning, which has since been fixed. The observed stellar 
continuum is ~ 1/4 the 2MASS flux levels for LkHa 233, reflecting the fact that our aperture does not include the whole extended source, 
possibly stellar variability and/or uncertainties in flux calibration. 

as well, but some authors have argued for a closer dis- 
tance; see ij3.4l below. Hernandez et al. also derived 



M* = 2.9 M Q , L* = 80 L Q , and age = 2.61 Myr, with 
Ay = 3.7, Ry — 5, but those values are likely not accu- 
rate, because they were derived without accounting for 
the fact that LkHa 233 is seen only in scattered light. 
Hence both the extinction to the star and its total lumi- 
nosity must be higher; see £13.21 

LkHa 233 sits at the center of an X-shaped bipolar 
nebula visible from optical to near -infrared wavelengths 
(jAspin et all 119851 : iLi et all [l99l . Near infrared po- 
larimetry shows a dark lane passing in front of the star, 
indicatin g the presenc e of a ro ughly edge-on cir cumstel- 
lar disk (|Perrin et ai1l2004fh ICorcoran fc Ravi (|1998af) 
discovered that LkHa 233 possesses an extended bipo- 
lar jet and counterjet visible in the optical [S II] 6716, 
6731 A lines. The blueshifted jet extends toward posi- 
tion angle (PA) 250°, perpendicular to the inferred cir- 
cumstellar disk and along the symmetry axis of the sur- 
rounding nebula. The redshifted counterjet, at 70°, was 
not detected within 0'.'6 of the star, which Corcoran et 
al. interpreted as oc cupation by the circum stellar disk. 
A followup study bv lMcGroartv et all (|2004D found that 
a chain of Herbig Haro objects extends almost 2 pc to 
either side of LkHa 233, located approximately along the 
jet axis but with considerable scatter. 

In the following section, §2, we present our observa- 
tions and outline the data reduction process. We then 
describe in §3 the observed properties of both the jet and 
the circumstellar disk, derive estimates for the mass flux 
in the jet, and re-examine the question of the distance 
to LkHa 233. In §4 we discuss the implications of these 
results, with special attention to comparison with the ob- 
served properties of jet around other Herbig Ae and T 
Tauri stars. Conclusions and prospects for future work 
are given in §5. 

2. OBSERVATIONS 

We observed LkHa 233 on 2005 Oct 13 UTC using 
the integral field spectrograph OSIRIS and laser guide 
star adaptive optics system (|Wizinowich et all 120061 : 
Ivan Dam et~al| |2006h on the W. M. Keck II telescope. 
The performance of the LkHa 233 itself served as the 



tip-tilt reference for the adaptive optics system. The 
skies were clear and seeing was average, 0'.'6 at K', with 
winds from 15-20 km hr -1 . 

OSIRIS is a lenslet-based integral field spectrograph 
which pr ovides a spectral re solution of 3800 from 1-2.5 
microns (|Larkin et al.ll2003f) . We selected its 50 mas 
pixel -1 spectrographic plate scale for our observations. 
While this scale does not Nyquist-sample the diffraction- 
limited PSF, it provides a larger field of view and better 
signal to noise per integration time for extended emission 
than the finer plate scales. We observed LkHa 233 us- 
ing the narrow-band HnS and KnS filters, covering the 
wavelength ranges 1.594-1.676 and 2.121-2.229 ^m, re- 
spectively. Total integration time on-source was 1800 s 
for HnS and 900 s for KnS. Between exposures on LkHa 
233, we obtained equally deep sky observations at a po- 
sition 20" away. 

We reduced the data using an early version of th e 
OSIRIS data reduction pipeline (jKrabbe et all l2004h . 
The data were sky-subtracted, corrected for offsets be- 
tween the 32 detector readout channels, spatially recti- 
fied to extract the individual spectra into a data cube, 
and corrected for cosmic rays and other blemishes. The 
A stars BD +9 5007 and HD 7789 served as telluric stan- 
dards and Point Response Function 2 (PRF) references. 
For BD +9 5007, observed immediately before LkHa 233, 
we measure a PRF FWHM of 61 ± 5 mas. We also used 
these stars to establish photometric zero points: we con- 
verted their 2MASS H and K s magnitudes into magni- 
tudes in the OSIRIS HnS and KnS filters via synthetic 

2 We follow the Spitzer convention of distinguishing be- 
tween the intrinsic optical Point Spread Function (PSF) 
and the Point Response Function (PRF), which is the 
PSF convolved with the detector's pixel response. See 
http : //ssc . spitzer . caltech . edu/postbcd/doc/PRF_vs_PSF . pdf 
Because our OSIRIS data used the 50 mas pixel scale, coarser than 
Nyquist for the Keck AO PSF at 1.6 fim, the PRF is dominated 
by the pixel response. By taking the quadrature difference of the 
observed 61 mas PRF FWHM and 50 mas pixel scale, we estimate 
that the intrinsic optical PSF had ~ 35 mas FWHM, approxi- 
mately the diffraction limit. Because of the difficulties inherent in 
trying to accurately measure Strehl ratios on undersampled data 
((Roberts et al. 200J), we do not attempt to directly measure the 
Strehl ratio achieved in these observations. 
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Fig. 3. — Continuum-subtracted velocity channel maps for [Fe II] 1.644 ^m emission from the jet. Each panel shows the intensity image 
for a given velocity (relative to the solar system barycenter) noted at the bottom of each panel. These are contiguous spectral channels 
from OSIRIS, and the field of view shown is 0'.'9 X 3'.'2. The blueshifted (bottom) jet is smooth and relatively featureless in both intensity 
and velocity, while the redshifted (top) jet breaks into three distinct clumps, with the central clump having a higher velocity than the other 
two. 



photometry, and then derived from those magnitudes a 
zero point for flux-calibrating our LkHa 233 data. 

These observations were made during the instrument's 
commissioning period, and show a variety of instrumen- 
tal artifacts, most notably ghost images at certain wave- 
lengths caused by crosstalk between lenslets due to a 
grating misalignment. These ghost images do not oc- 
cur in the wavelength ranges of interest to us, and so 
we simply ignore them. The grating has been realigned 
subsequent to these observations and more recent data 
do not suffer from these ghosts. 

OSIRIS also includes an imaging camera whose field of 
view is offset by 19'.'4 from the spectrograph field. Our 
sky observations were obtained such that LkHa 233 fell 
on the imager field of view (FOV) while we were obtain- 
ing spectrograph skies, and vice versa. In this manner 
we obtained i?, J, H, and K band observations of LkHa 
233, with total exposure times of 300, 360, 360, and 315 
s, respectively. These observations were reduced via the 
usual steps of flat fielding, sky subtraction, image regis- 
tration, and summation. 

During this observing run, we also observed the Her- 
big Ae/Be stars LkHa 198, V376 Cas, and Parsamian 
21, with identical instrument configuration and expo- 
sure times. No outflows were visible around the first two 
sources. Parsamian 21 has a faint outflow extending to 
the north, which we postpone discussing for a subsequent 
paper. 

3. RESULTS 
3.1. The Jet 

A well-collimated bipolar emission-line jet extends sev- 
eral arcseconds on either side of LkHa 233 (Figure [1]) , 
with the blueshifted jet at a position angle of 249°, and 
the redshifted jet at 69°. We also see a roughly edge-on 
circumstellar disk and bipolar nebula extending perpen- 
dicular to the jet. The jet has excavated a cavity in this 
envelope, the limb-brightened edges of which form the 
well-known "X" centered on LkHa 233. We detect the 
jet in the [Fe II] 1.600 and 1.644 /im lines, and possi- 
bly H 2 2.122 (im (Figure [2|). Faint hydrogen Br7 2.16 
/im emission is also detected along the course of the jet, 
but does not show velocity shifts and is spread over a 
wider spatial extent than the [Fell] emission; this ap- 
pears to just be Br7 emission from the star, scattered 
by dust in the circumstellar envelope. Of the lines ob- 
served here, the 1.644 /im line is the brightest, and we 



choose to concentrate our analysis on it. In Figure [3] we 
show continuum-subtracted images in a range of velocity 
channels centered on 1.644 /im, and we show in Figure 
|4] a position-velocity diagram, an summed continuum- 
subtracted jet image, and a plot of the FWHM. 

The jet is narrow, tightly collimated, and knotty. Its 
FWHM increases slowly from 100 AU near its base to 
~ 200 AU at the edge of our FOV, 1300 AU from the 
star, giving an opening angle of 9° (Figure [4j. Three 
bright knots are apparent in the redshifted jet, at approx- 
imate distances of 0.3, 0.8, and 1.3" from the star, cor- 
responding to 260, 700, and 1150 AU in the plane of the 
sky. We will refer to these as knots C, B, and A, respec- 
tively. The blueshifted jet is noticeably less clumpy, but 
brightens considerably near its base, reaching its peak 
intensity at O'.'l (88 AU) from the star. There is an an- 
ticorrelation between jet width and intensity (Figure |4j 
right), with the bright knots having lower FWHM than 
the fainter inter-knot emission. This behavior has also 
been observed for the outflows from RW Aur (Woitas et 
al. 2002) and HH 23 (Ray et al. 1996). 

We were unable to find any reported measurements 
of LkHa 233's radial velocity in the literature 3 , nor are 
there stellar lines suitable for measuring it in our data 
(since Br~/ is in emission). Hence we compute velocities 
relative to the solar system barycenter, without correct- 
ing for the radial motion of the star. In this reference 
frame, the red- and blueshifted jets have average veloc- 
ities of +108 ± 3 and —89 ± 3 km s _1 , respectively. If 
the red and blue jets have the same speed relative to the 
star but opposite directions, then the stellar barycentric 
radial velocity appears to be ~ — 10 km s^ 1 . 

A significant difference between the two directions is 
that the blueshifted jet has a relatively constant velocity 
along its length, while the redshifted jet shows variation 
in its velocity. Knot B, at 0'.'8, has a higher velocity 
than the other two clumps by about 30 km s" 1 (Figure 
IU left). This velocity variation so close to the jet's origin 
indicates that the amount of acceleration imparted to the 
jet must fluctuate with time, and must do so in a way 
that is not symmetric between the jet and counterjet. 

If we assume the jet is perpendicular to the disk, which 
has an estimated inclination of i = 65° ± 5 (see i )3.2l 

3 We note that Corcoran and Ray (1997,1998) state they mea- 
sured the systemic velocity using optical Na D lines, but do not 
anywhere report the actual measured velocity 



Bipolar Jet from LkHa 233 



5 




Fig. 4. — Left: Continuum-subtracted position velocity diagram for [Fe II] 1.644 fim, using a O'.'l wide "virtual slit" centered on the 
jet. Center: Continuum-subtracted image of the jet, integrated in velocity. The upper and lower portions of this image are summed over 
slightly different wavelength ranges centered on the red- and blueshifted jets, each with bandwidth 0.6 nm. The vertical dotted lines show 
the boundary of the "virtual slit" used to produce the left-hand figure. Right: Measurements of jet FWHM versus distance from the star. 
This plot was constructed by fitting each row in the continuum-subt racte d image with a Gaussian. The FWHM appears narrower at the 
position of the knots, and broader where the emission is fainter (see 93.lt . 



below), then the radial velocities we observe are only 
40 ± 9 percent of the total jet velocity. The implied total 
jet velocity, 200-300 km s" 1 , is co nsistent with outflow s 
seen around other Herbig Ae stars (jMundt fe R ay 1994). 
Based on this, we estimate the proper motion of the jet 
should be 0'.'05 ± 0.01 per year. Hence the bright knots 
in the jet were launched from the star within the last 10 
to 20 years, and their motion should be apparent over 
short timescales. 

We do not see much evide nce for the low ve l ocity c om- 
ponent (LVC) reported by ICorcoran fc Ravi (1998a) at 
about -25 km/s. This may be due to the fact that the 
[Fe II] 1.644 /im line traces higher densities than the [S II] 
6716, 6731 A lines, so in the NIR we are only sensitive 
to dense gas in the core of the jet. 

The emission lines of [Fc II] can be used to measure 
several important physical quantities, s uch as temper- 
ature , electron density, and extinction (jHartigan et al.l 
200 1;. Most such measurements require the compari- 
son of flux ratios between different lines of [Fe II] , or 
between [Fe II] and other species in the optical. The 
limited wavelength coverage of our data restricts their 
diagnostic capabilities, as only the 1.644 and 1.600 /.tm 
lines are present in our wavelength range. The flux ra- 
tio of those two lines does provide a measure of th e 
electron density in the shocked gas (jNisini et al.l 12002). 
We find that the electron density in LkHa 233's jet is 
10 4 - io 4 - 5 cm and decreases at larger distances 
from the star; see Figure [5] This decrease at larger dis- 
tances has also bee n seen in many jets f rom T Tauri stars 
dNisini et all l2002t iPodio et alj 120061: iBeck et all [20071 ; 
lHartigan fe Morsell2007ft and is believed to show the slow 
recombi nation of hydrogen after the initial shock ion- 
ization (|Bacciotti fe Eisloffell [l999h . However, as [Fe II] 
emission traces denser regions of the gas, our measured 
electron densities are not necessarily indicative of LkHa 
233's outflow as a whole. Future observations which 
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Fig. 5. — The ratio of the 1.644 fim and 1.600 fim lines pro- 
vides a diagnostic for electron density n e . The solid and dashed 
lines show the expected ratio as a function of n e for two differ- 
ent temperatures, from Nisini et al. (2002). The measured line 
ratios for the three clumps in the redshifted jet of LkHa 233 are 
indicated; the length of the crosses shows the uncertainties. The 
electron density decreases along the jet as one moves away from 
the star. The values of ~ 10 4 - 10*' 5 cm 3 are comparable to 
electron densities in sev eral T Tauri HH jets me asured from the 
same [Fe II] line ratio bv lNisini et~aTl ((2003. 12005T) ■ 

also include the optical [S II] doublet are needed to pro- 
vide complementary measures for gas at lower densities 
and/or different distances from the shock fronts. 

3.2. The Disk 

In addition to showing the jet, these observations also 
clearly reveal the circumstellar disk around LkHa 233, 
in the form of a bright nebula crossed by a dark lane, 
the characteristic signature of a roughly edge-on circum- 
stellar disk (Figure [TJ right). The presence of this disk 
was originally sus p e cted based on optical p olarimet ry 
(jAspin et al.lll985fl . ICorcoran fe Ravf (|1998al ) were un- 
able to trace the redshifted optical [S II] emission any 
closer than 0'.'6 to the star, and inferred the presence of 



Perrin and Graham 



0.5 I 



Observed 



• 





Model Image 




10" 




SED 












0.5 


















10 










e 






^ Model 


0.0 






10 13 


/ 






I 1 














> 






Input Star 








10 14 






0.5 


















:[ 


\ \ 








10 15 







-0.5 



0.0 
Arcsec 



0.5 



-0.5 0.0 0.5 

Arcsec 



1.1 1.0 10.0 100.0 1000.0 
M^rn) 



Fig. 6. — Left: Observed appearance of the disk around LkHc* 233. This 1.6 /im image is the median of the continuum from the 
OSIRIS Hn3 datacube. Center: A 1.6 fim disk model computed with MCFOST, convolved with a measured OSIRIS Hn3 PRF obtained 
from observations of a calibration star. The model reproduces the central brightening, sharp V-shaped southern disk face, and less distinct 
northern disk face. Right: The observed SED for LkHa 233 (diamonds), plus the computed SED from our model (solid line). The SED is 
fairly well fit for wavelengths < 50 fim; the deficit at longer wavelengths arises because the model presented here includes only the disk and 
not the surrounding extended nebula, which contains cooler dust at much larger distances from the star (but still within the large IRAS 
60 and 100 /mi beams). The dashed line is the input stellar spectrum, and the solid line shows the resulting spectrum after reprocessing 
by the disk. See §3.2 in the text for a discussion of this model and its limitations. 



a disk with this radius that occulted the jet. We confirm 
here that interpretation: the circumstellar disk extends 
from (X'5 below to 0'.'6 above the locus of peak intensity. 
In contrast to the [S II] line, in [Fe II] 1.644 light 
the jet can be traced all the way in to the star on both 
sides (FigureE]). There is at most a two-pixel-wide (O'/l) 
region at the base of the red jet where the emission is 
somewhat fainter, possibly due to obscuration. This in- 
dicates that the circumstellar dust is not optically thick 
at these wavelengths, in contrast to the optical. The ob- 
served disk morphology is consistent with that seen in 
visible light in HST WFPC2 observations (Stapelfeldt et 
al., in prep). 

The continuum intensity peak is well-resolved in these 
data, with a FWHM of tf/135 versus Of! 06 for a PRF 
reference star. Therefore we do not observe the star di- 
rectly, but instead see it in disk-scattered light. The peak 
intensity is located at the apex of the southwest face of 
the disk, which is the brighter of the two faces. This is 
consistent with the southwest side of the system being 
inclined toward us, as indicated by the blueshifted jet. 
By comparison with a scries of circumstellar disk mod- 
els computed us ing the MCFOST Monte Carlo radiative 
transfer code of iPinte et al.l (|2006f ). we determined that 
LkHa 233's disk is inclined approximately 65±5 degrees 4 . 
A numerical model of a circumstellar disk at this incli- 
nation extending from 0.6 to 500 AU reproduces fairly 
well both the observed appearance and SED of LkHa 
233 (Figure [6|). We caution that some parameters, such 
as the inner disk radius, are not well constrained by these 
data. Hence this is only one possible disk model which 
fits the current data, not necessarily a unique best-fit 
solution. 

Previous estimates of the stellar properties of LkHa 
233 itself were obtained assuming only regular inter- 
stellar extinction, typically arou nd Ay = 2 — 4, with 
reddening Ry bet ween 3 and 5 ([Hernandez et al.1 12004 
iManoi et aLll2006[ ). However, our results here confirm 

4 We measure inclinations of the disk midplane relative to the 
plane of the sky, such that pole-on disks have i = 0° and edge-on 
disks have i = 90° 



that LkHa 233 is seen entirely in scattered light at opti- 
cal and near-IR wavelengths. The true Ay is much higher 
than previously estimated, so the observed SED contains 
only a small fraction of the total luminosity. LkHa 233 
must therefore be substantially more luminous than pre- 
viously thought. With the stellar temperature fixed at 
T» = 8600if (appropriate for an A5 star), we were able to 
simultaneously match both the observed disk morphol- 
ogy and SED by using a stellar radius of i?» = 9R© . This 
implies a total stellar luminosity ~ 400L e — about an 
order of magnitude higher t han t he 20-80 L©previously 
found by iHernandez et al.l (|2004h . These revised stel- 
lar properties place LkHa 2 33 near the stellar birthline 
for intermediate mass stars (Palla fe Stahler| [l990; Palla 
2005), and therefore LkHa 233 may in fact be an very 
young star with mass ~ 4M Q . 

The key parameters of our model are as follows: T* = 
8600 K, i?* = 9 R Q , r out = 500 AU, and grain sizes from 
0.03 — 10 ^jm with power law index -3.5. However, we em- 
phasize that the model presented here is not a rigorous 
best fit to the data, merely one plausible fit, and there 
remain large uncertainties in both disk and stellar prop- 
erties. In particular, the MCFOST Monte Carlo code 
used does not yet account for any puffe d up inner rim 
to the disk (|Dullemond fe D ominik 2004) nor docs it in- 
clude the extended bipolar envelope also present aroun d 
LkHa 233 in addition to the disk (|Perrin et al.1 12004T) . 
The lack of the bipolar envelope component results in the 
present model having an unrealistically large scale height 
(30 AU at 100 AU) in order to fit the data well. More de- 
tailed modeling, quantitatively fit to observations across 
a wide range of wavelengths, would help clarify the situ- 
ation, but is beyond the scope of this paper. 

3.3. Mass Flux 

The observed intensity in the [Fe II] lines allows us 
to estimate the mass flux. Currently the physical pa- 
rameters of the jet, such as the temperature, ionization, 
and Fe depletion onto grains, are not well known so our 
estimate will require making several assumptions; fu- 
ture observa t ions w ill enable a more refined calculation. 
iPodio et al.l (|2006f ) presented two separate methods for 
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estimating mass flux based on [Fe II] emission, based on 
either (A) the jet's apparent size and density, or (B) the 
observed emission line luminosity. 

The first method relies on computing the mass outflow 
rate across a given plane perpendicular to the jet, based 
on the observed jet cross section, velocity, and density: 

M = fj,mHnHTTr^vj 

where /i is the mean molecular weight, taken to be 1.24, 
tuh is the mass of hydrogen, tih is the number den- 
sity of hydrogen, and vj is the jet velocity. In the case 
of LkHa 233, we do not know njj, but we do know 
that n e ~ 5 x 10 4 cm -3 , from which we can estimate 
n,H- Typical ionization fractions for HH objects are 
x e = 0.03 — 0.3, as determined fr om both shock mod- 
els and observational dia gnostics (lHartigan et ah) 119941 : 
iBacciotti fe Eisloffel|[l999tlLavallev-Fouquet et aUboOOj j. 
so if we assume x e = 0.1, then Hh — 5 x 10 5 cm -3 . 
Using a rough estimate of rj = 90 AU, we compute 
M = 1.0 x 10 -6 M Q year -1 . This method relies on the 
assumption that the jet can be treated as a uniformly 
filled cylinder with density measured by our diagnostic, 
and thus may be better considered an upper limit on the 
true M. 

The second method relies on the observed luminosity 
L\Feii} being proportional to the total number of emit- 
ting atoms. Unlike the first method, this second ap- 
proach does take into account spatial variations in the 
jet density, but is subject to uncertainties in the absolute 
flux calibration. In this method we compute M using 



M = iim H (n H V)v t / k 



(n H V)=L [FeII] I hvAif, 



Fe+Fe 
Fe IT 



where V is the volume filled by the emitting gas, Vt is the 
tangential velocity of the jet. It is the tangential length 
of the jet volume element, ^jpt and 4? are the ioniza- 
tion fraction and atomic abundance relative to hydro- 
gen, Ai is the Einstein A radiative rate coefficient, and 
h and v are as usual the Planck constant and frequency 
of the radiation. Based on our estimated disk inclina- 
tion i — 65° ± 5 and observed v ra diai = 100 km s , we 
compute v t ~ 240 ± 40 km s -1 , and we take l t to be the 
tangential length of each OSIRIS pixel, 44 AU at 880 pc. 
Several assumptions are needed in order to transform the 
observed line luminosity £[_F e /j] into the density- volume 
product njjV: without accurate knowledge of the tem- 
perature we have to estimate the ionization fraction and 
upper state fraction for iron. We base our estimates on 
conditions o b served in ty pical HH objects, fo r insta nce by 
iPodio et all (|2006f) and lHartigan fe Morsel (|2007l) . who 
measured temperatures of 8000 — 25000 K based on opti- 
cal line ratios. Taking those temperatures plus our mea- 
sured n e = 10 4 — 10 4,5 . we u sed a 16-level model of the 
Fe+ ion (Pcse nti et al.l 120031 . similar to that described 
by) to estimate the fraction of ions in the upper energy 
level is — 0.01 ± 0.005. We further assumed that Fe 
is entirely ionized (Ep£- — 1), and that iron has a solar 
abundance but is 90% depleted onto solid grains (based 
on the 87% depletion reported for the HH 23 jet by Podio 
et al. 2006). 

Our spatially resolved IFU data enable us to measure 
the mass flux rate as it varies along the jet. We first 



collapse the 2d continuum-subtracted [Fe II] image per- 
pendicular to the jet to get a total £[Fe//] for each row, 
and then apply the above algorithm to compute M. The 
computed mass flux has a median value of 1.2 x 1O -7 M0 
year -1 , and is fairly constant along most of the jet, vary- 
ing between 1.0 — 1.5 x 10 -7 M Q year -1 . The one excep- 
tion is the middle knot in the redshifted jet, knot B, 
where the flux rises to 2.6 x 10 -7 M Q year . We caution 
that these values are approximate at best, given the lack 
of knowledge of the temperature and ionization struc- 
ture, and furthermore represent lower limits, because all 
of the outflow is not necessarily shocked and luminous, 
and because we have not corrected for extinction. 

The two methods give results which differ by an or- 
der of magnitude, with the estimate based on method A 
higher than that of method B. This discrepancy is not en- 
tirely unexpected, since the filling factor for [Fe II] emis- 
sion is less than unity. Podio et al. (2006) found masses 
computed with method A were often 3-10 times higher 
than those computed based on emission line fluxes, and 
they considered the latter method to provide a better es- 
timate of the mass flux transported by the atomic com- 
ponent of the outflow. 

The values for M for LkHa 233 are comparable to the 
high end of the range of mass fluxes seen around T Tauri 
stars. For instance, if we adopt the lower of the two esti- 
mates as the more reliable, the flux in LkHa 233's jet is 
roughly comparable to the mass flux computed for RW 
Aur by Woitas et al. (2002), ~ 1O -7 M year -1 , though 
about an order of magnitude higher than the values mea- 
sured for HH 34 and HH 111 by Podio et al. (2006). 

Outflow rates from young stars are proportional to 
the accretion rates onto those stars (|Cabrit et al.1 119901 : 
lHartigan et al. l 1995: E dwards et al.ll2006f). across a mass 
range spanning from b rown dwar fs ( Moha ntv et al.|[2"005f) 
to Herbig Ae/Be stars (jCorcoran fe Ravll998bl ). Outflow 
mass flux es are typically a few percent of the accreting 
mass flux (|Hartigan et al.lll995l ) As far as we know, there 
are currently no measurements of the accretion rate onto 
LkHa 233, but based on our estimated outflow mass flux, 
the accretion rate of onto it should be ~ 1O 5 M year -1 . 

3.4. The Distance to LkHa 233, Revisited 

As discussed above, the fact that LkHa 233 is seen en- 
tirely in scattered light implies that it must be substan- 
tially more luminous than has previously been thought, 
assuming it is indeed at a distance of 880 pc. It is worth 
considering whether it might instead be closer. 

Establishing accura te distances to Herbig Ae stars is 
often problematic (e.g. iStine fe O'Nea l 1998; Merm et al.l 
12004 IPerrin et al.ll2006h . The 880 pc distance t o LkHa 
233 was originally derived by ICalvet fe Cohenl (|1978f) . 
and the majority of subsequent authors have adopted 
their estimate. Calvet & Cohen derived that value from 
the distance modulus to the nearby star HD 213976, a 
bright Bl star located 6.'5 away on the edge of the dark 
cloud in which LkHa 233 is embedded. They also es- 
timated photometric distances to two other members of 
this small cluster, the K stars LkHa 231 and 232, finding 
d = 600 — 700 pc subject to the assumption that those 
stars are both overluminous by 2 mags, in V due to their 
youth. 

Is it possible that HD 213976 is in fact not at the 
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same distance as LkHa 233? Aspin (1985) noted that 
HD 213976's optical colors imply its Ay is only 0.42, 
significantly less than the typical Ay = 2.5 per kpc. 
This implies both that we look toward it along a galac- 
tic sightline with lower than average extinction, and that 
HD 213976 must be located on the near side of the dark 
cloud. Hence LkHa 233 should be at least as far away as 
HD 213976, if not fa rther. 

Some authors (e.g. lBowev et al.ll2003[ ) have instead de- 
scribed LkHa 233 as being a member of condensation A 
of the Lac OBI molecular cloud, w hich is only 418 ± 15 
pc distant (jHernandez et al.ll2005[ K While LkHa 233's 
position on the sky does place it on the apparent edge of 
Lac OBI, we consider it unlikely to actually be a mem- 
ber of that association, for several reasons. First, LkHa 
233's proper motion is very different from that of Lac 
OBI : (pi cos b,fib) = (—9,21) mas/yr (jDucourant et al.1 
2005) vers us (m; cos b, Hb) = (—2 .3 ± 0.1, —3.4 ± 0.1) for 
Lac OBI; (|de Zeeuw et al.lll999D . Even more problem- 
atic is_tfiejfe£t_that 'Lac OBI is believed to be 16 Myr 
old (|de Zeeuw et al.|[l~999D . Because protoplanetary disks 
dissipate on time scales of 3-10 Myr, the presence of a 
massive disk around LkHa 233 is inconsistent with its 
being that old. For these reasons we prefer to retain the 
traditional 880 pc distance for LkHa 233. 

Because the system's inclination can be constrained 
based on the disk's appearance, the proper motion of the 
clumps in the jet can potentially provide a direct measure 
of the distance. For a disk inclination i = 65° ± 5 and 
d = 880 pc, the jet proper motion should be 0'.'06 ± 0.01 
year -1 , while d = 420 pc would imply a proper motion of 
0'.'12 ± 0.03 year -1 . However, this will be a challenging 
measurement given uncertainty in the inclination, the 
possibility of jet precession, variations in the internal 
working surfaces traced by the clumps, etc. 

4. DISCUSSION 

The outflow from LkHa 233 is in general quite sim- 
ilar to those seen around T Tauri stars. The velocity 
structure and collimation of the jet, the morphology of 
the knots, and the presence of asymmetries between the 
jet and counterjet all resemble the properties of outflows 
seen around lower-mass stars. 

For instance, we see that for both the red and blue 
directions, the higher velocity channels arc more tightly 
collimated than the lower velocity channels. Such in- 
creased collimation at higher ve locities has also bee n 
seen in the outflow from DG Tau (|Bacciotti et al.ll2002f ). 
This is in agreement with predictions of MHD models 
that inner streamlines should dominate the emission (e.g. 
iDougados et al.ll2004h . 

It is striking how similar LkHa 233 appears to the 
T Tauri star HH 30, with its well-known edge -on disk 
and perpendicular outflow (|Burrows et al.|[l996t ). LkHa 
233's disk is asymmetric in brigh tness between the left 
and r ight sides, as is HH 30's disk (W atson fc Stapelfeldtl 
120041) . The outflows from both stars are knotty, and even 
the characteristic periods for knot production are similar, 
2.5 years for HH 30 and w 5 years for LkHa 233. These 
knots appear to be due to internal shocks within the jet, 
rather than the jet shocking into the ambient medium. 
The properties of the knots in LkHa 233's jet are indeed 
consistent with internal working surfaces caused by vari- 
ations in jet outflow velocity. In particular, the bright- 



est knot is also that with the highest velocity, just what 
would be expected for more rapid material plowing into 
a slower portion of the jet. 

It is curious that, at our angular resolution and S/N, 
only the redshifted jet appears knotty. Similar asym- 
metries have b een seen in outflo ws from about half of 
T Tauri stars (|Hirth et al.lll994h . On the other hand, 
at least s ome jets are very sym metric, such as HH 212 
in Orion (jZinnecker et al.l lT998). which displays a very 
regular series of knots in nearly perfect pairs on opposite 
sides of the star. The symmetry of the HH 212 knots and 
the regularity of knot spacing in many systems (including 
LkHa 233) suggest that the production of knots is an in- 
herent part of the jet launch process, possibly due to disk 
instabilities or recurrent stellar activity cycles. However, 
the commonness of asymmetric systems indicates that 
any inherent regularity in the process can easily be dis- 
rupted or masked, by collisions with surrounding ambient 
material, pressure gradie nts, or asymmetries in magnetic 
fields very near the star (jWassell et al.l 12006). But then 
what causes those asymmetries in magnetic field or pres- 
sure? The origin of these various symmetries and asym- 
metries remains mysterious, and poses a key challenge to 
theories of jet acceleration and collimation. 

A particularly useful comparison object is HD 163296, 
which has perhaps the best-studied outflow from 
any Herbig Ae star (|Devine et al.l 120001 : IWassell et al.1 
l lOOg). HD 163296 an A2Ve star at 122±^ pc 
(|van den Ancker et al.1 Il998h . launches a bipolar colli- 
mated jet, HH 409, which also is broken into an asym- 
metric series of knots. The HH 409 jet displays knots 
on both the jet and counterjet, but on the two sides the 
knots differ in both and distance from the star. Based on 
their high spatial resolution observations of HD 163296, 
Wassell et al. concluded that the asymmetries must orig- 
inate near or within the jet launch region itself, rather 
than being the result of interactions with the ambient 
medium. Other authors have reached the same conclu- 
sion regarding the asymmetries observed in outflow s from 
T Tauri stars (jWoitas et alj|2002t iPvo et al.ll2006t ). Sim- 
ilarly for LkHa 233 , given that strong asymmetries are 
present in the outflow within the innermost few hundred 
AU, it seems unlikely that those asymmetries could be 
due to the ambient medium. While there is indeed sub- 
stantial dust and gas near the star (in the observed disk 
and nebula), we know the outflow from LkHa 233 has 
been ongoing for thousands of years based on the length 
of the observed outflow (at least 1.3 pc; McGroarty et 
al. 2004). Therefore any primordial material along the 
jet axis was long ago swept away, out of the inner few 
arcseconds where we observe the asymmetries and knots. 
Explaining the difference in appearance between the red- 
and blue-shifted jets thus requires the presence of an 
asymmetric geometry within the inner jet launch region. 
This asymmetry may take the form of pressure gradi- 
ents, different magnetic field strengths or differing field 
geometries between the two directions. 

The mass flux in the jet of LkHa 233 is about two 
orders of magnitude greater than that of HD 163296. The 
measured electron density n e for LkHa 233 is also about 
a factor of ten higher than that reported for HD 163296, 
but the value for HD 163296 was measured using ratios of 
optical [S II] transitions which trace lower density parts 
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of the post-shock gas (|Wassell et al.l [20061 ). Hence it is 
hard to say how meaningful this difference in n e is, and 
the true difference in accretion rate may be somewhat 
less than this estimate. 

Outflow mass fluxes are believed to be proportional to 
accretion rate, a nd accretion rates are generally larger 
for younger stars (|Hartmann et alj [l998) . Therefore, to- 
gether with its higher extinction, greater infrared excess, 
and larger amount of nebulosity, the high mass flux of 
LkHa 233 suggests that it is at a younger evolutionary 
state than HD f 63296. This is consistent with our revised 
stellar parameters which place LkHa 233 near the stel- 
lar birthline for intermediate mass stars. The outflow 
velocity is expected to be of order the escape velocity, 
v esc = y/GMn/R*. For our revised stellar properties (4 
M Q , 9 Rq), the escape velocity is v esc w 300 km s _1 , 
which compares favorably with the estimated total jet 
velocity. 

There is a growing body of evidence that outflows from 
high mass stars are generally less collimated than those 
of low mass stars. Given the tightly collimated nature 
of LkHa 233's outflow, the transition to poorly colli- 
mated jets appears to occur at a higher mass. However, 
the opening angle of the bipolar cav ity in the envelope 
around LkHa 233 (jPerrin et al.ll2004D is 30°, much wider 
than the jet collimation angle. This, plus the scatter in 
HH objects seen by McGroarty et al. (2004), argues ei- 
ther for a change in outflow direction with time (such as 
from precession), for a much lower degree of collimation 
during earlier stages of the outflow, or for the presence 
of a wide-angle wind around the inner jet, unseen in the 
[Fe II] emission which traces relatively dense gas. 

5. CONCLUSIONS 

We have used Keck Observatory's laser guide star 
adaptive optics system and integral field spectrograph 
OSIRIS to observe the bipolar outflow of the Herbig Ae 
star LkHa 233. Based on these observations, we believe 
that: 

1. LkHa 233 has a narrow, collimated bipolar outflow, 
with the blueshifted jet at a position angle of 249° 
and the redshifted counterjet at 69°. 

2. At the present angular resolution and sensitivity, 
the blueshifted and redshifted jets have similar 
opening angles but very different appearances, with 
the red-shifted jet broken up into several knots 
while the blue-shifted jet is smooth. The knots 
in the red-shifted jet are consistent with shocked 
inner working surfaces within the jet. 

3. Based on [Fe II] emission intensity, we estimate the 
mass flux in this jet to be (1.2 ± 0.3) x 10 -7 M© 
year -1 . This is toward the top end of mass fluxes 
observed around other T Tauri and Herbig Ae 
stars, and suggests that LkHa 233 is a very young 
source with strong ongoing accretion. 

4. The observed properties of LkHa 233's outflow - 
its kinematics, collimation, asymmetric and knotty 
appearance, electron density, and estimated mass 
flux - are all largely consistent with the ob- 
served properties of jets from T Tauri stars. The 
transition to the less-collimated outflows observed 



around massive YSOs must therefore occur at some 
still higher mass threshhold. 

5. The jet extends perpendicular to a flattened cir- 
cumstellar disk, the innermost region of the X- 
shaped nebula around LkHa 233. 

6. LkHa 233 is seen only in scattered light from 
this disk; the star itself remains obscured at near- 
infrared wavelengths. Because of this, previous es- 
timates of the extinction to the star are low, caus- 
ing the total stellar luminosity to be underesti- 
mated in prior studies. 

7. Revised stellar parameters, based on fitting disk 
models to the SED, indicate that LkHa 233 may 
be a 4 M Q star near the stellar birthline, but there 
remain large uncertainties in its true stellar prop- 
erties. 

The main result of this study is confirmation that the 
outflow from LkHa 233 is quite similar to typical outflows 
from T Tauri stars, as has been suggest ed based on prior 
observations at lower spatial re solution ([Corcoran fe Ravi 
Il998al iMcGroartv et al.ll2004l ). This supports the hy- 
pothesis that the same physical mechanism is responsible 
for creating stellar outflows, across the entire mass range 
from T Tauri to Herbig Ae stars. It is worth noting that 
circumstellar disks around Herbig Ae stars are also sim- 
ilar in broad outline to those around T Tauri stars (e.g. 
Gradv et aD2005t IPerrin et alJl2ot)5lMillan-Gabet et all 
2007|) In contrast, there are suggestions that circumstel- 
lar disks around the Herbig Be and more massive YSOs 
differ from those around lower masses. In particular, 
there is evidence for a transition from magnetic accretion 
(for T Tauri and Herbig Ae stars) to direct disk accre- 
tion onto the star (for Be st ars and beyond) (|Vink et al.l 
120021 iMottram et al.l l2007f ). Neither this transition in 
accretion mechanisms nor the transition from collimated 
to wide-angle outflows are particularly well understood. 
But the fact that both transitions occur at approximately 
the same stellar mass may be an important clue that 
these processes are connected. This may indicate either 
that the properties of outflows depend primarily on the 
properties of the circumstellar disks and how they accrete 
onto the star, or that both outflow and disk properties 
depend in turn on some third parameter (stellar mag- 
netic fields?) which undergoes a transition at the Ae/Be 
boundary. 

Further insight into these matters will depend on ob- 
taining increasingly detailed spatially resolved spectra of 
outflows, such as those provided by integral field spec- 
troscopy. We have shown here the measured electron 
density using the 1.644/1.600 fim ratio, and have also 
presented a rough estimate of the mass flux using the 
1.644 /im line. But much more could be done with a 
more complete spectrum: Ratios of other lines (most 
notably 1.644/1.533) provide alternate measures of elec- 
tron density, and comparison with o ptical [Fe II] lines 
can measure the temperature as well (|Nisini et al.l 120051 : 
iPodio et al.ll2006f K The 1.644, 1.320, and 1.257 /zm lines 
all originate in the same upper state, and thus their ra- 
tios directly constrain the extinction Ay . This is a par- 
ticularly exciting prospect for the redshifted jet in LkHa 
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233: it is located partially behind the inclined circumstcl- 
lar disk (occulting it in the visible but not in the near 
infrared) , and thus it can be used to measure the extinc- 
tion through the disk. Future observations that include 
the complete optical and near-IR spectrum of [Fe II] and 
other shock tracers should enable a complete suite of di- 
agnostics for temperature, density, ionization, refractory 
element depletion, and more. Detailed modeling of LkHa 
233 thus holds great promise for understanding outflow 
mechanisms around intermediate-mass stars. 
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dous efforts of the OSIRIS instrument team and Keck 
Observatory staff, whose members are too numerous to 
list here in full. Special thanks in particular go to 
James Larkin, Shelley Wright, Mike McElwain, Randy 
Campbell, and Al Conrad. Conor Laver assisted with 
the OSIRIS data reduction pipeline. Christophe Pinte 
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